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Abstract 

In-medium nucleon-nucleon scattering cross sections are explored by comparing results of 
quantum molecular dynamics simulations to data on stopping and on elliptic and directed flow 
in intermediate-energy heavy- ion collisions. The comparison points to in- medium cross sections 
which are suppressed at low energies but not at higher energies. Positive correlations are found 
between the degree of stopping and the magnitudes of elliptic and directed flows. 
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One of the main goals of research in the area of heavy ion collisions (HICs) at intermediate 
energies has been the determination of bulk properties of nuclear matter, such as the nuclear 
equation of state (EoS). While a considerable progress |5| including that in 

subthreshold kaon production which provides complementary information on the EOs[(| 
has been reached in determining the EoS at supranormal densities, relying on reaction data, 
the uncertainties are still very large. To access the EoS, it is necessary to describe reaction 



observables 
reaction theory 



such as those quantifying the collective motion of nuclear matter, within 
1 1 if ] . The transport models employed in the description 



IS, Q, y, y, 
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of central reactions have included the quantum molecular dynamics approaches in its QMD 



nnn 



Q, H3| and ImQMD y,M,H6| (with 
Uehling-Uhlenbeck (BUU) approaches 



m for Improved) variants, as well as the Boltzmann- 



17 



191 ] . The two main ingredients of the nuclear 



transport are the nucleonic mean fields and nucleon-nucleon binary scattering cross sections 
(NNCS). The employed cross sections affect virtually any observable from a central reactions 
and constraining those cross sections is essential for reducing the EoS uncertainties 9J. Also 
the in-medium cross-sections are of interest for their own sake, as they underly the viscosity 
and other nuclear transport coefficients j20j ]. Best for constraining the cross sections are 



the stopping observables 
However, the correlation 



21 



22] that reflect the cross sections in the most direct manner. 



21] with flow observables, such as used for determining the EoS, 



is also of interest. In t 
ImQMD model Q, ljj 



lis paper, for constraining the cross sections, we shall employ the 



16|]. 



Principally, within theory the mean field and cross sections of transport calculations 
should be linked to the same microscopic nucleon-nucleon interactions. In practice, this has 
turned out to be overly ambitious resulting in the phenomenological strategies in transport 
of varying independently the mean field, related to EoS, and the NNCS. Microscopically, far 
fewer calculations have been done of the in-medium cross-sections than of the EoS and of the 
mean fields. Specifically, the zero-tem per ature in- medium elastic NNCS (ENNCS) have been 



studied wit 
InRefs. 2J 



lin the relativistic 
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24 



251 ] and non-relativistic Brueckner approaches [26|, [2 



251 ] . employing the Brueckner relativistic approach, it was found that, compared 
to the free ENNCS, the in- medium ENNCS were suppressed at low relative momenta and 
less suppressed at higher relative momenta and even enhanced slightly depending on the 



medium density. The in-medium ENNCS have also 



path Green's function (CTPGF) approach [28, 



29 



3een studied by using the closed time 



30 



31] employing, in particular, the 
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QHD-I and QHD-II effective Lagrangians, with the mean field and in-medium EN1N 



derived self-consistently for the same effective interaction. The latter studies 



29 



30 



CS 



3l| 



produced in-medium correction factors which were different for the cross sections of like, 
°nriDD an d unlike, cr* F , nucleons. Otherwise, the zero-temperature ENNCS from the Dirac- 

3l| have exhibited similar qualitative features but 



nn,pp 

Brueckner 



3, 



25| and CTPGF 



29 



30 



the results have differed in the details of their dependance on density and energy. The 



temperature dependence of the in-medium ENNCS has been investigated in 



31] within 



the CTPG 
Elsewhere 



approach and a general increase in the ENNCS with temperature was found. 



32j , though, difficulties have been pointed out in defining the cross sections within 



a nuclear medium. Those difficulties might relegate the cross sections utilized in transport 
to strictly phenomenological quantities. In the comparisons of nuclear reaction simulations 
to data, indeed, most often phenomenological cross section parameterizations have been 
employed, as e.g. represented by the formula 



a 



NN 



WPoVnn 



free 



(1) 



where i] = 0.2 
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341 ] . Also, the in-medium ENNCS scaled by the effective mass, 



^nn/^nn = (m*(p,p)/fn) 2 has also been employed in comparisons [35| to data from heavy- 
ion collisions (HIC). The latter scaling presumes that, for given relative momenta, the matrix 
elements of interaction are not changed between the free space and medium. 

The purpose of this work is to draw conclusions on the in-medium ENNCS using recent 



data on stopping and elliptic and directed flows, obtained with a good centra! 



from collisions of Au + Au and other symmetric or near-symmetric systems 



21 



ity selection, 
39, ba. For 



those data, a high degree of correlation was found between the degree of stopping and the 



strength of collective flow. 



the ImQMD model 



HQ 



n our investigations we rely on the recent version ImQMD05 of 



16(]. 



Within the ImQMD05 model, the mean fields acting on nucleon wavepackets are derived 
from an energy functional where the potential energy U includes the full Skyrme potential 
energy with just the spin-orbit term omitted: 



U — Up + U m d + C/coui 



(2) 



Here, ?7coui ^ s the Coulomb energy, while the nuclear contributions can be represented in 
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local form with U p>Tad = J d 3 r u Pjia d and 



UL H . y r i i/sur \2 , i/sur,iso r-r-* / \ii 



m - ^ 4 - p7+1 
p 2 po 7 + 1 Po ' 2 Po Po 

+ (V + 5p 7+1 + Cp 8 / 3 )5 2 + ^4^' ( 3 ) 

Po 

where 5 = (p n — p p )/(p n + Pp), P = Pn + P P and p n and p p are the neutron and proton 
densities, respectively. The energy associated with the mean-field momentum dependence 
may be represented as 

u md = ^- Yl 7^6 / dVd 3 p 2 /iY 1 (p 1 )/iV 2 (P2)1.57 [ln(l. + 5. x l(T 4 (Ap) 2 )] 2 , 

(4) 

where /jy are nucleon Wigner functions, Ap = |px — p 2 |, the energy is in MeV and momenta 
are in MeV/c; the resulting interaction between wavepackets is such as in Ref. 38j. The 



coefficients utilized in ([3]) can be transcribed, see [161 ] . onto those usually specified for the 
Skyrme interactions. In this work, the SkP and SLy7 Skyrme interactions are employed. 
Both of those interactions give rise to an incompressibilty of K ~ 200 MeV and produce an 
EOS consistent with the features of collective flow in HICs from Fermi to relativistic energies. 
However, the symmetry energies associated with those two interactions are different. It 
should be stressed that the expression for the energy ([3]) is more complete here than in the 
preceding Ref. 

Within this paper, three different phenomenological forms of in-medium ENNCS are 
utilized. The first set cr NN are the cross sections given by Eq. (|TJ) with p = 0.2 and with 



the free cross sections described in terms of the parameterization of Ref. 39 
set cr™ are cross sections calculated within the CTPGF approach of Ref. 



The second 



301 ] . following 



the QHD-II Lagrangian. For calculational convenience, within ImQMD the formula ([T]) is, 
though, employed with r\ made dependent on density and energy and its dependence fitted 
to the CTPGF results. For reference, in Fig. 1 we present both the cross sections, o*^ 
and ann,pp in panels (a) and (b), and the suppression parameters, rjnj and r]^ pp ^ in panels 
(c) and (d), as a function of cm. energy a/s at different densities. One can see that 
changes little with density and is nearly the same as in free space. On the other hand, the 
cross section <Tn^ PP tends to be suppressed at lower energies, y/s < 2.05 GeV and enhanced 
at higher. Differences in the features of the two cross sections are associated with the 
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differences between the T = and the T = 1 channels and, in particular, presence of a low- 
energy resonance in the T = 1 channel and effects of p exchange. The third set constitutes 
and ad hoc parameterization, inspired by the CTPGF results, aiming at the description of the 
excitation function for elliptic flow in the midrapidity region of |y C m/2/cm am | < 0.1 in Au + Au 
collisions In that parameterization, the common cross-section modification parameter, 
^( 3 ) = rjnj = f]np , depends on the beam energy for reaction, -Ebeam, in the following manner: 
r]^ = 0.2 for £ bcam < 150AMeV, r)® = for 150AMeV < £ bcam < 200 AMeV, 77® = -0.2 
for 200 AMeV < E hcam < 400 AMeV, and r/ 3 ) = -0.4 for 400 AMeV < £ beam . In this paper, 
we actually confine ourselves to the HIC energy range of ^beam < 400 AMeV since, on one 
hand, our model is nonrelativistic and, on the other, higher energies require a consideration 
of the inelastic cross section that we are not prepared to carry out at this moment. 

First, we investigate the impact of in- medium NNCS onto the model predictions for ellip- 



tic and directed flows in Au + Au collisions. 

m 



have been determined experimentally [21 



Excitation functions for both flow observables 
The top panels in Fig. 2 show the ex- 



citation function of the elliptic flow parameter t>2 for Z < 2 particles in the midrapi dity 
region of \y C m/ycm m \ < 0-1- The v% values have been obtained using the rotated frame [36 ]. 
The bottom panels in Fig. 2 show the excitation function of maximal scaled directed flow 
in Au + Au collisions at b ~ 5fm (&/6 max ~ 0.38). The scaled flow is defined with [2lJ 
P«L = PsdirACS where u h ™ m = 7 c b ^ am ^ am and p xdil = £ sgn(y) Z u x / £ Z. The sums in 
the directed flow extend over all particles with charge number Z < 10. The data and calcu- 
lations are represented in Fig. 2 with closed and open symbols, respectively. Lines generally 
merely guide the eye, except for the dashed lines in the bottom panels. In these cases, the 
dashed lines illustrate the change in the directed flow data after applying the correction for 
reaction-plane fluctuation. The calculations in the left panels of Fig. 2 have been carried 
out with the SkP and SLy7 mean fields and a*^ cross sections. The differences between 
the calculated V2 values in the panel (a) of Fig. 2 are limited and can be due to the differ- 
ence in the symmetry energy between the two mean-field models. Notably, there are some 
differences between the data sets in the figure too. One more significant difference which 
develops between the calculations for different mean fields in Fig. 2 is in p^ iT in panel (c) 
at .Ebeam > 1 50 AMeV. The calculations in the right panels of Fig. 2 are for the SkP mean 
field in combination with the three different cross-section models. The sensitivity of flow 
to the cross sections is similar to the sensitivity to the mean-field models. It is limited 
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on absolute scale but gets more enhanced at higher energies. The deviations between the 
models are similar to the deviations of models from data. From the presented mean-field 
and cross-section combinations, the SkP mean-field in combination with a*^ cross sections 
describe the flow data best. 

We next turn to the impact of NNCS on nuclear stopping. Recently, in the exper- 
iment 21|, the ratio of the rapidity variance in the transverse direction to the rapidity 



variance in the longitudinal direction, vartl, has been used a measure of the nuclear stop- 
ping. The longitudinal rapidity for the ratio is defined in the standard manner in the cm 
system. The transverse rapidity is defined by replacing the longitudinal direction in the 
definition with a random transverse direction. Figure 3 shows a variety of results pertinent 
to vartl. The results are for charged particles with Z = 1 — 6 from central collisions of 
symmetric or near-symmetric systems, with the contributions of different particles weighted 
with Z. The panels (a) and (b) in Fig. 3 show calculated distributions in longitudinal and 
transverse rapidities in 400 AMeV Au + Au collisions at b/b maiX < 0.15. The corresponding 
vartl values for different calculations are quoted in those panels. The panel (c) compares 
chosen calculated Au + Au excitation functions for vartl to data. Finally, the panel (d) 
compares the calculated dependence of vartl on system charge to data at 400 AMeV. 

From the different panels of Fig. 3, the panel (a) tests the sensitivity of the calculated 
rapidity distributions to a utilized mean field. It is apparent that that sensitivity is quite 
weak. Quantitatively, when switching from the SkP to Sly7 mean field, vartl increases just 
by 0.015 i.e. relative 2%. The panel (b) next tests the sensitivity of the calculated rapidity 
distributions to a utilized cross section. The sensitivity is much greater here, with vartl 
changing by 0.22 or 32% when switching from the a*^ to a*^ cross section. Thus, there 
is a good chance to restrict the in-medium cross sections using measured vartl but less 
chance to restrict the mean field. What is important in Fig. 3(b) is that the larger the 
cross section in the given energy region the more similar are the transverse and longitudinal 
distributions and the larger the corresponding vartl value. The panels (c) and (d) compare 
to data the results obtained with the SkP mean field and two of the in-medium cross sections, 
cr*^) and a*^ . At lower energies a semiquantitative agreement with the Au + Au data is 
found and there the cross sections coincide. However, at higher energies the Au + Au data 
favor a*W NNCS. The <7*« NNCS appears excessively reduced and this also concerns a*^ > 
given the 400 AMeV vartl value of 0.805 in the panel (b). Looking next at the system 
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charge dependence of vartl in the panel (d), it appears that an optimal in-medium cross at 
400 AMeV could be a bit lower than a* 1 - 3 ). Of the three cross sections, nonetheless, a*^ 
compares best with the data. 

Experimentally, the values of vartl maximize at i?beam ~ 400 AMeV and so do the values 
of pj^ir- At the same time, the elliptic flow values minimize in this energy region 36]. If the 
values of vartl at different energies for the Au + Au system are considered simultaneously 
with the values of maximal scaled directed flow, a positive relatively narrow correlation is 
found between the observables 2JJ. Importantly, that correlation is between observables 



determined, on one hand, in the most central b/b max < 0.15 collisions and, on the other, 
the semicentral b/b max ~ 0.38 collisions. This suggests, from the experimental side, that 
both observables probe bulk nuclear properties which come into play at specific densities 
and excitations reached at a particular beam energy. Between the vartl and v 2 values, again 
determined at different centralities, a negative correlation is expected. Interestingly, when 
other symmetric or near-symmetric systems are considered, and specifically Ni + Ni, Ru + 
Ru and Xe + Cs, the resulting correlations appear to extend the correlation lines found for 
Au + Au, see Fig. 4, with the left panel (a) showing the p^^-vartl correlation and the right 
panel (b) the V2-vartl correlation. The system systematics appears further to underscore, 
from the experimental side, that bulk nuclear properties get tested by the observables. 

It is obviously important to find out whether the ImQMD05 model produces correlations 
comparable to those measured and whether the calculated and measured correlations agree 
quantitatively. Fig.4(a)and (b) show the correlation of maximal p^ ir and v% with vartl, 
respectively. All of these results are for the reaction systems Ni + Ni, Ru + Ru, Xe + 
Cs,Au + Au at E beam = 250, 400AMeV. While the data are represented with crosses 
and stars in Fig. 4, the calculations are represented with planar-figure symbols. For each 
of the represented cross sections, i.e. either a*^ represented by open symbols, or a*^ 
represented by closed symbols, in combination with the SkP mean field, the calculations 
produce approximate correlation lines in the planes of flow vs vartl, with positive slope in 
the case of p^ ir and negative in the case of t>2. One of the calculated correlations, V2-vartl 
for a*^\ is somewhat broader than the others. To aid the eye, the correlation areas for a*^ 
are additionally marked with shadowing. Within the reaction simulations, the correlations 
between stopping and flow may be primarily associated with the collisions governed by cross 
sections. The greater the momentum transfer in collisions the higher the stopping. However, 
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the occurring local randomization due to collisions also enhances hydrodynamic behavior 
and, thus, enhances the magnitude of flow observables. From the calculated correlations 
in Fig. 4, only the two obtained using o~*( 3 ' reasonably agree with data. Notably, with the 
stopping and flow observables simultaneously growing in magnitude with collision number, 
the correlation lines associated with the different cross sections could have actually coincided. 
However, vartl and flow observables depend on interactions differently. The vartl observable 
exhibits a strong dependence on cross sections and little on mean field, while flow observables 
depend to a comparable extent on cross sections and mean field. As a result, the correlation 
lines can end up being different for the different cross sections. 

In summary, within the ImQMD05 model we have investigated the impact of in-medium 
NNCS and mean fields on stopping and flow in HICs at intermediate energies. The vartl 
observable, quantifying the stopping, depends strongly on NNCS and little on mean-field 
details. On the other hand, the v 2 and p^ ir flow observables depend comparably on un- 
certainties in the NNCS and mean fields. Within the upper range of the energies we have 
investigated, both the common NNCS parameterization represented by a*^ and the CTPGF 
cross section represented by a*^ appeared excessively reduced in the medium to produce 



observed vartl values 21] . Also the magnitudes of flow observables calculated for those 
cross sections turned out to be low compared to the experiment. The agreement between 
the data and calculations gets much improved when employing an ad hoc cross section pa- 
rameterization guided by CTPGF and represented by a*^ 3 \ For that parameterization the 
nucleon-nucleon cross sections in the medium are suppressed at low energies and enhanced 
at higher energies, i.e. the extracted in-medium nucleon-nucleon cross sections depend on 
both density and energy. In this respect it is also in agreement with the prediction of 241.125]. 



The results clearly do not represent the last word on in-medium cross sections as medium- 
variation of the cross sections in the parameterization has been clearly oversimplified. From 
the side of observables, those that we examined do not strongly differentiate between the 
like and unlike nucleon cross sections. 
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Fig.l Energy dependence of in-medium cross-sections and cr^pp, in panels (a) and 

(b), and of the suppression parameters r/np and ^nn, PP , in panels (c) and (d), at selected 
densities. The cross sections have been obtained within the CTPGF approach with 
QHD-II effective lagrangian. 

Fig. 2 Flow excitation functions in Au + Au collisions. The top panels display midrapidity 
elliptic flow. The bottom panels display maximal scaled directed flow. The data are 
represented by solid symbols and calculations by open symbols. Lines generally merely 
guide the eye, except for the dashed lines in the bottom panels. In these cases, the 
dashed lines illustrate the change in the directed flow data after applying the correction 
for reaction-plane fluctuation. The left panels in the figure illustrate the sensitivity of 
calculations to the employed mean field. The right panels illustrate the sensitivity of 
calculations to in-medium cross sections. 

Fig. 3 Comparison of characteristics of central colliding systems in the longitudinal and 
transverse rapidities. Panels (a) and (b) display the calculated rapidity distributions 
of particles from central (b/b max < 0.15) 400 AMeV Au + Au collisions and illustrate, 
respectively, the sensitivity of calculations to the mean field and to the cross sections. 
The rapidities are normalized to the beam cm rapidity. Panels (c) and (d) display 
the calculated and measured dependencies of the variance ratio vartl, respectively, on 
energy for the Au + Au collisions and on net system charge at collision energy 400 
AMeV. The data are from the FOPI Collaboration 2l|. The calculations have been 



done for the SkP mean field combined with either o"*^ 1 -* or a*^ cross sections. 
Fig. 4 Correlation between vartl and flow observables: p^j lr in the panel (a) and v 2 in (b). 



Both the data 
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